Hybrid Electron Spin Resonance and Whispering Gallery Mode Resonance 

Spectroscopy of Fe 3+ in Sapphire 



Karim Benmessai^'E] Warrick G. Farr, 1 Daniel L. Creedon, 1 Yarema 
Reshitnyk, 2 Jean-Michel Le Floch, 1 Timothy Duty, 3 and Michael E. Tobar 1 

1 ARC Centre of Excellence for Engineered Quantum Systems, 
University of Western Australia, 35 Stirling Highway, Crawley WA 6009, Australia 
2 ARC Centre of Excellence for Engineered Quantum Systems, 

University of Queensland, St Lucia QLD 4067, Australia 
3 ARC Centre of Excellence for Engineered Quantum Systems, 
University of New South Wales, Sydney NSW 2052, Australia 
(Dated: March 15, 2013) 

The development of a new era of quantum devices requires an understanding of how paramagnetic 
dopants or impurity spins behave in crystal hosts. Here, we describe a new spectroscopic technique 
which uses traditional Electron Spin Resonance (ESR) combined with the measurement of a large 
population of electromagnetic Whispering Gallery (WG) modes. This allows the characterization 
of the physical parameters of paramagnetic impurity ions in the crystal at low temperatures. We 
present measurements of two ultra-high purity sapphires cooled to 20 mK in temperature, and 
determine the concentration of Fe 3+ ions and their frequency sensitivity to a DC magnetic field. 
Our method is different to ESR in that it is possible to track the resonant frequency of the ion from 
zero applied magnetic field to any arbitrary value, allowing excellent measurement precision. This 
high precision reveals anisotropic behaviour of the Zeeman splitting, which has not been previously 
reported. In both crystals, each Zeeman component demonstrates a different ^-factor. 



Since the development of the first Whispering Gallery 
maser oscillator based on Fe 3+ ions in sapphire, the 
system has been the focus of research at several 
institutions^— However, the physical parameters of these 
paramagnetic ions in sapphire have never been charac- 
terised with sufficient accuracy. Several experiments have 
attempted to measure the spin-spin relaxation time, or 
the number of ions involved in the maser process, but 
have given results which contradict the values given by 
the manufacturer of the crystals, or those found in the 
literature.— ~— This difference is essentially due to the 
method used to interact with the ions in sapphire, namely 
the excitation of high-Q Whispering Gallery modes at 
frequencies coincident with the ESR frequency of the ion. 
The spatial distribution of these modes acts as filter, se- 
lecting only the population of ions which fall inside its 
field pattern. Furthermore, the inhomogeneous broad- 
ening of the ion's ESR bandwidth makes the characteri- 
zation of their physical parameters difficult. Until now, 
values in the literature have not been confirmed by an- 
other method^*^ The present work describes a new tech- 
nique for measuring these parameters, which in principle 
is valid at any temperature. In our work, crystals are 
cooled to near 20 mK in a dilution refrigerator, and all 
WG modes within a few GHz of the ESR frequency of 
Fe 3+ in sapphire (12.04 GHz) are identified. Then, a 
magnetic field is applied parallel to the cylindrical z-axis 
of the crystal, and swept in value. At each value of mag- 
netic field, the reflection and transmission coefficients of 
the WG mode population are measured. Characterisa- 
tion of the frequency and quality factor of the modes 
allows the ESR bandwidth and the concentration of the 
ions to be determined accurately enough to confirm the 
work published previously. This step is very important 



for such a system, because it is an excellent candidate for 
use in a host of quantum devices, with sapphire typically 
containing many species of residual paramagnetic spins 
compared to the current engineered quantum systems. 
Due to the inhomogeneous broadening effect, different 
spin packets arc dispersed to different frequencies. It has 
been shown that these may act independently, with a 
coupling between the spin packets governed by the hyper- 
fine interaction with the aluminium in sapphire lattice^. 
These 'self interactions have also been observed exhibit 
non-linear behaviour. Combining this with the high Q- 
factors of the WG modes, sapphire is an attractive ma- 
terial for future quantum devices J^r— It is therefore im- 
portant to fully characterize the physical parameters of 
any ion defect centers in the crystals. 



I. ION DESCRIPTION 

The fundamental energy levels of the Fe 3+ ion in sap- 
phire are defined by the zero- field splitting parameter 6 5 
in the spin Hamiltoniam 12 ' 13 ! 18 ' 19 Three levels exist at 
zero applied magnetic field, with the dependence of these 
levels on a DC magnetic field described by the following 
Hamiltonian: 
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where g « 2 is the Lande g-factor, /is = 9.274 x 10 
J T _1 is the Bohr magneton, B is the DC magnetic field 
strength and S is the electron spin angular momentum of 
the ion. The Hamiltonian parameters are summarized in 
Table UJ These parameters were measured at 4.2K dur- 
ing the 1960s and vary from crystal to crystal due to 
different spatial distribution and concentration of ions 
in the lattice, and due to crystal dislocations and varia- 
tions in the symmetry of the crystalline field. Here, our 
calculations use the parameters given by Kornicnko and 
Prokhorov as they give results close to those measured 
in our crystals. The Hamiltonian parameters have never 
been calculated for extremely low temperatures, near 20 
mK. Nevertheless, we assume in our calculations that 
they do not change significantly from their values near 
4 K, especially the g-factor, transition frequencies, and 
sensitivity to DC magnetic field. 





Symmons 


Kornicnko 


D 


1719.2 ± 1 


1838.5 ± 0.6 


\a\ 


229.4 ± 1 


253.5 ± 1.3 


a-F 


351.5 ± 1 


362.7 ± 2 



TABLE I. Hamiltonian parameters for Fe 3+ in sapphire, 
given by Kornienko & Prokhorov and Symmons & Bogle at 
4.2 K in units of 1(T 4 cm" 1 . 

Zccman splitting occurs in the presence of a DC mag- 
netic field, with each energy level split into two sub-levels 
as shown in Fig. [I] Each sub-level has a sensitivity dif- 
ferent to the others, and crossing of sub-levels occurs 
at large values of magnetic field. Allowed transitions 
between sub- levels follow the selection rule Am = ±1. 
At zero magnetic field the ions occupying the first en- 
ergy level will split in to two distinct components; one 
which we denote ESRi ow corresponding to the transition 
|+l/2) — > |+3/2), where the frequency or energy differ- 
ence decreases when B increases, and the other which 
we denote ESRhigh- This corresponds to the transition 
| — 1/2) — > | — 3/2) where the frequency or energy differ- 
ence increases with B. 

In this work, we examine the transition between the 
lowest energy levels. The choice of these transitions i 
justified by the setup condition (Section [IT]), and because 
nonlinear behaviour has previously been observed at their 
frequencies. The energy difference is expressed for the 
two transitions as a function of the DC magnetic field 
strength B (expressed in Tesla) as follows: 

AE l±1/2) _^ ±3/2) (B) = ±B - D + ^(a - F) 

+ 7^(95 + IBD + (a - F)) 2 + 80a 2 (2) 
o 

and the corresponding transition frequencies are ex- 
pressed as 

^high/low = — 7^^-^+1/2) -+±3/2) ■ (3) 
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FIG. 1. Zeeman splitting of the ESR energy levels of Fe 3+ in 
sapphire in the presence of a DC magnetic field. The presence 
of the sapphire cubic field, represented by the a term in the 
Hamiltonian, creates a small admixing between the | ±5/2) 
and | + l/2) states (e ~ 0.03)^. 

The predicted evolution of the transition frequency be- 
tween the states is shown in Fig. [2] 
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FIG. 2. Evolution of the frequency of the transition between 
Zeeman levels with B. 



II. EXPERIMENTAL SETUP 

The crystals used are a cylindrical HEMEX-grade sap- 
phire from Crystal Systems (USA), 5 cm diameter x 3 
cm height, and containing 1-3 parts per million (ppm) of 
Fe 2+ and Fe 3+ impurity ionsi22r— In the present work, 
we present results using two crystals, denoted CI and 
C2, in which the Fc 3+ concentration is measured to be 10 
and 100 part per billion (ppb), respectively.^ The crys- 
tals are cut such that the c-axis of the sapphire lattice 
is parallel to the cylindrical z-axis. To attain the high- 
est Q-factor at low temperature, the crystal is ordinarily 
mounted in the center of a cylindrical silver plated cop- 
per cavity. Energy is coupled to the resonator through 
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two straight antennas, placed 180° from each other, en- 
tering through the top of the cavity and oriented par- 
allel to z-axis. The energy is stored in the crystal in 
high azimuthal order Whispering Gallery modes where 
the H field is predominantly oriented perpendicular to 
the z-axis. The crystal was cooled using a dilution re- 
frigerator (DR), with a cooling power of 1 W at the 4 
K stage, and ~300 //W at 100 mK. A superconducting 
5 Tcsla magnet was mounted to the 4 K stage of the 
DR such that samples connected to the mixing cham- 
ber plate (~20 mK) could be placed in the bore of the 
magnet. A sapphire was mounted in the center of the 
magnet coil, with the walls of a radiation shield used 
to confine the energy of the WG modes rather than the 
usual cavity, which was too large to fit into the bore of 
the magnet. As a result, the Q-factors of the modes 
were strongly degraded and reached values around 10 6 
instead of 10 8 to 10 9 which is typical with a high quality 
cavity. Measurements of the WG modes were performed 
using a network analyzer with an excitation power 
of -40 dBm at the input of the probe antennae. This 
power was low enough to allow measurements to be made 
without perturbing the temperature of the system. At 
higher power levels, the temperature of the sapphire rises 
quickly to ^250 mK with an associated frequency shift. 
Using low levels of applied power also allowed the ob- 
servation of absorption effects due to the presence of 
paramagnetic impurity ions, and allowed complex ther- 
mal bistability effects to be avoided, which have been 
observed in such a system in the pasti^ 3 . When the input 
power level is extremely low, all the energy inside the 
crystal is absorbed, making the observation of the ESR 
variations with the application of a magnetic field easy 
to do. For higher power measurements, the temperature 
was controlled at 400 mK. 

As described earlier, the transitions of interest were 
|±l/2) — >• |±3/2). The choice of characterising these 
transitions only was due to several experimental limita- 
tions. The network analyzer used for the experiment was 
limited to measurements below 20 GHz, and in addition, 
the microwave transmission lines installed in the dilu- 
tion fridge presented very high loss above 18 GHz due 
to their connector type. The cables were connected via 
fecdthroughs between each temperature stage of the DR, 
with attenuators connected at each feedthrough to allow 
full thcrmalisation of the cables and reduce temperature 
gradients in the system. The crystals were characterized 
up to a magnetic field of 600 G, which corresponds to a 
frequency interval of 10-14 GHz for the ESR of Fe 3+ . At 
higher magnetic field, more careful measurement must 
be made because of the complexity of the Zeeman spec- 
tra. The energy levels of different ions cross each other 
and it becomes difficult to identify them clearly with- 
out following the 19 GHz line upwards from G. If we 
take in account ESR of other ions present in the sapphire 
such Cr 3+ , Mn 4+ and V +2 r^2£ the system quickly be- 
comes very complex to study, even if the concentration 
of these ions is very low. The magnetic field value was 



set manually through a magnet controller with a ramp 
time of about one minute between field values. After the 
set value was reached and the persistent current switch 
of the magnet closed, measurements were only performed 
after the temperature of the system had stabilized at the 
base temperature of ~20 mK. 
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FIG. 3. Schematic of the experimental setup. 

III. MEASUREMENT METHODOLOGY 

The characterisation of the ions in our crystals is 
based on measurements of the AC susceptibility of the 
WG modes around the Fe 3+ ESR frequency when a DC 
magnetic field is applied parallel to the c-axis. The 
S-parameters of the WG modes between 10 and 14 GHz 
are measured, which act as filters sharply defined in fre- 
quency that allow the frequency of the ESR to be deter- 
mined as the applied magnetic field is changed. When 
the frequency of the ESR is far from a WG mode, no 
change is observed. When the ESR frequency is coinci- 
dent with a WG mode frequency, the S-parameters will 
decrease in amplitude, the frequency of the mode shifts, 
and the Q-factor of the mode is significantly degraded. 
This process is fully described in Section llll Al and related 
to the AC susceptibility of the ions. The lineshapes of 
the WG mode resonances are modelled with a Fano reso- 
nance fit. The equation of the fit is shown in Fig. 31 and 
allows the mode frequency and Q to be computed. From 
the results obtained, the AC susceptibility is calculated, 
and the physical parameters of the ions are determined. 

A. AC magnetic susceptibility 

The magnetic susceptibility of the ions in sapphire is 
described by the following general equation: 



x(f) = x» + *x"M 



(4) 



We assume here the ESR is homogeneously broadened, 
and the resonance lineshapes are Lorentzian. As will be 
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FIG. 4. Example traces of a WG mode showing the transmis- 
sion S-parameter for zero applied DC magnetic field, and for 
a 32 Gauss applied field. The WG modes are degenerated and 
show generaly a double resonance due to cavity imperfections 
and alignements 
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FIG. 5. Real and imaginary parts of the magnetic suscepti- 
bility as a function of magnetic field strength for one of the 
WG mode doublet at 12.390 GHz. 



demonstrated in Section IIIIE1 the ions are in fact inho- 
mogencously broadened at zero applied magnetic field. 
Here, we assume them to behave homogeneously when 
the Zeeman components of the ESR are split with the 
application of a magnetic field and are far away from 
each other in frequency. The analytical equation of the 
AC magnetic susceptibility is given by Kramers-Kronig 
relationships: 



1 + (2ttt 2 ) 2 {v - v tJ ) 2 + -r 1 r 2 ( 7 i7) 2 ( 7 2 



(5) 
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1 + (2ttt 2 ) 2 (v - v v} ) 2 + -TiT 2 (jH) 2 a. 



where, n and t 2 are the spin-lattice and spin-spin relax- 
ation times of the transition ij, respectively, afj = 2 for 
the transitions at 12 GHz, the term ^T]T 2 (^iB) 2 a 2 ^ <C 1, 
and the DC magnetic susceptibility is xo- The real and 
imaginary parts of the magnetic susceptibility are related 
to the frequency shift and Q-factor of the modes as fol- 
lows: 



( v WG (B = Q)-v WG {B) 1 , 
v B 2 



VX 



(6) 



I QY G (B) QY G (B = 0) 



where, vwg arL d are the frequency and loaded Q- 

factor of the WG mode. As an example of the fit (Fig. [5]) 
shows the real and imaginary parts of the magnetic sus- 
ceptibility of one mode at vwg=^2.390 GHz. At low 
magnetic fields, the Fe 3+ ESR is centred around a sin- 
gle frequency. When the magnetic field increases, the 
ESR splits into two Zeeman components with two reso- 
nant frequencies: vi ow , which shifts to the lower frequen- 
cies, and vuigh to which shifts higher in frequency toward 



12.390 GHz - the frequency of the WG mode considered 
here. The frequency of the mode, vwg, will increase 
as Vhigh approaches it, until it reaches a maximum and 
decreases to return to the initial value, at which point 
Vhigh = vwg- When the magnetic field is further in- 
creased, Vhigh continues to increase and the frequency of 
the WG mode decreases to a minimum before returning 
to its initial value. The Q-factor of the mode stays con- 
stant while Vhigh vwg i but as it becomes close, the 
losses due to the presence of the ions increases and the 
Q- factor of the modes decreases. When Vhigh 3> vy/G-, 
the Q-factor then increases to reach it's initial value. The 
characterisation of all WG modes between 10 and 14 GHz 
allowed the measurement of the frequency dependence of 
each Zeeman component. Figure |6] shows the evolution 
of the ESR frequency with magnetic field strength. The 
center frequency is defined as that where the real part of 
the susceptibility of the mode crosses zero, or when the 
imaginary part reaches it's minimum. 



ESR frequency and sensitivity to DC magnetic 
field 



As described earlier, the AC susceptibility allows the 
frequency of the Fe 3+ ESR to be determined when sub- 
jected to an applied DC magnetic field. By following the 
shift of the Zeeman components, it's possible to obtain 
their frequency sensitivity Shigh/iow to the magnetic field, 
as shown by Fig. |6] The theoretical value is about 2.8 
MHz/Gauss, assuming isotropic Zeeman splitting. How- 
ever, the results of the fit show that this sensitivity is 
also dependent on the concentration of Fc 3+ ions in the 
sapphire. Indeed, for both crystals S'high is found to be 
higher than S\ ovn and the sensitivity of the sapphire (C2) 
is higher than for the second crystal, CI. The theoretical 
prediction of 2.8 MHz/Gauss is obtained from the first 
order approximation of Eqns. [5] and [3] by using parame- 
ters shown in Table |U This fact highlights the novelty of 
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this work because the characterisation of the paramag- 
netic ion in the past was performed with the traditional 
ESR method where a microwave signal is applied to the 
crystal in a microwave cavity at a fixed frequency near 
the ESR of the ion, and then a DC magnetic field is 
applied to the cold system. The measurement of the fre- 
quency shift allows the calculation of the real part of the 
ac susceptibility and then the estimation of the Hamil- 
tonian parameters and sensitivities. In our case, we map 
out the path of the ESR frequency from zero field to 
the desired value of the magnetic field, giving a more ac- 
curate estimation of the Zeeman components frequency 
sensitivities. Our results reveal the anisotropic nature of 
these components, and we later give an estimation of the 
p-factor based on these results. 
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FIG. 6. Frequency sensitivity to DC magnetic field for two 
sapphires at 20 mK. Only the points for a magnetic field 
higher than 20 Gauss are considered because of the complex- 
ity of the ions behavior (cf. section HlIE[) . 

The frequencies as a function of DC magnetic field are 
given in Fig. [51 The fit shows that the zero-field ESR fre- 
quency is different for the upper and lower Zeeman com- 
ponents, which is consistent across both crystals. The 
crystal C2 contains has an ion concentration 10 times 
that of CI, and it has been shown in a previous publica- 
tion how the greater concentration of ions revealed very 
complex behaviours. In the present case, the suscepti- 
bility obtained for low magnetic field (B < 20 Gauss) 
does not follow a Lorentzian shape and is not symmetric 
around the mode frequency, making the fit invalid and 
estimation of the correct resonance frequency difficult. 
As a first estimation, the central frequency is situated in 
the range 12.028-12.050 GHz. We expect that the two 
plus and minus spin packets do not perfectly overlap at 
zero magnetic field and form a large ESR. The different 
sensitivities are show in Fig. [7] where the measured val- 
ues are compared to the parameters in in Table |TJ The 



sensitivities are presented as a function of the ESR fre- 
quencies at zero field. It is clear that for both crystals the 
sensitivity of the ESRhigh is higher than the sensitivity 
of the ESRiow Also, the sensitivity of the crystal with 
a higher ion concentration is greater than for the lower 
concentration crystal. 
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FIG. 7. ESR Zeeman component sensitivities as a function of 
frequency. 



C. Determination of g-factor 

The measurement of g is well known, and different 
techniques have been employed in the past for that pur- 
pose. A one-electron quantum cyclotron measurement is 
the most accurate, with a precision exceeding one part 
per trillion, measuring g/2 = 1.00115965218085(76)322 
The determination of g has been achieved in numerous 
materials and techniques, such GaAs, CdTe, InP, and 
graphite^ - — In AI2O3 crystals, electron spin resonance 
measurements have allowed the calculation of this factor 
for different ions in the microwave range. These mea- 
surements are summarised in Table [ill Other measure- 
ments of energy spectra in the optical and millimetre 
wave regimes have allowed the determination of this fac- 
tor with great accuracy. For the Fe 3+ ion in sapphire, 
the g-factor determination is a poorly represented mea- 
surement in the literature.— 

TABLE II. Measured (/-factor for a variety of impurity ions 
and host materials. 



Material 


Ion 


(/-factor 


Reference 


A1 2 :! 


Fe y + 


2.0026±0.0005 
2.0034±0.0003 


V2 

13 


A1 2 3 


Mo 3+ 


1.98±0.01 (_L) 
1.968±0.001 (||) 


2k 


A1 2 3 


Cr 3 + 
(ruby) 


1.9867±0.0006 (_L) 
1.9840±0.0006 (||) 


Stt 


Ti0 2 


Fe y + 


2.000±0.005 


Jil 


Ti0 2 


V 4+ 


ff*=1.914 

ffy=1.912 

^=1.956 


38 


Emerald 




1.97±0.01 (_L) 
1.973±0.002 (||) 


39 



Kornicnko & Prokhorov, and Symmons & Bogle have 
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FIG. 8. Lande g-factor vs. frequency as measured in this work 
in two sapphires, CI and C2, as well as the values measured 
by other authors. The g— factor is extracted from the fit of the 
WG modes frequency sensitivity to DC magnetic field as one 
parameter. 



estimated this parameter and claimed that g is indepen- 
dent of temperature, while they found the other Hamil- 
tonian parameters to be temperature dependent. This is 
in contradiction to other materials and ions where the g- 
factor is ordinarily temperature dcpcndcnt ) 30 i 33 thus we 
expect that the earlier measurements were not achieved 
with enough accuracy to reveal temperature induced vari- 
ations. In addition, g is expected to be isotropic for Fe 3+ 
in sapphire in the direction perpendicular to the crystal 
c-axis. In the present work, the g-factor is obtained from 
one parameter linear fit of the data points shown in Fig. 
[6l Figure [8] shows the g- factor as a function of the reso- 
nance frequency. For the two crystals, g is smaller for the 
upper Zeeman component than for the lower component. 
The values are far from the those already estimated in 
the past. However the mean value of g is close to that ob- 
tained by Bogle and Prokhorov. In addition, it appears 
that the g— factor is dependent on the ion concentration. 
Where the concentration is higher, the g— factor is closer 
to 2.003 for both Zeeman components. The g-factor is 
also affected by an important parameter of the crystal, 
which is the alignment of the cylindrical z-axis with the 
crystal lattice c-axis. HEMEX-grade crystals show a tol- 
erance of 1° that affects the ^-factor by Ag = 0.00015, 
which can be added to the fit tolerance quadratically. 



D. Magnetic losses due to ions 

The WG modes of the resonator are characterised by 
the Q-factor given by Eqn. 042^ The loaded Q-factor 



of the mode (QY G ) is related to the unloaded Q-factor 
(Qo), which is a function of the external losses (Q e ) and 
the losses due to the presence of magnetic ions (Q m )- 
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1 



(7) 



The magnetic quality factor Q m is related to the mag- 
netic losses in the resonator, and can be determined by 
the decreasing value of Q^ G with the application of a 
magnetic field. When the frequency of the ESR does not 
match the mode frequency, then the magnetic losses are 
negligible (Q^ 1 — > 0). However, when the ESR frequency 
is tuned equal to the WG mode frequency, this magnetic 
loss mechanism becomes important and the Q-factor of 
the mode can be loaded significantly. 



1 

Qm 



l 



:) 



(8) 



The magnetic Q is related to the magnetic loss tangent 
tan((5 m ) as follows 



1 

Qrr, 



77tan(r5 m ) 



(9) 



The tan(<5 m ), or magnetic loss tangent, is shown in Fig. [S] 
as a function of frequency at a temperature of 20 mK. The 
losses due to the presence of the Fe 3+ ions are similar for 
both Zeeman components, and stay relatively constant 
around 2 x 10~ 5 . 
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FIG. 9. Magnetic losses due to Fe 3+ as a function of fre- 
quency. For clarity reasons, only the most relevant results of 
crystal CI are presented. 



E. Inhomogeneous broadening effect 

At zero applied DC magnetic field, the lineshape 
broadening of the Fe 3+ ESR is not well known. Previous 
work has shown that the ions are inhomogencously broad- 
ened, which has resulted in the observation of a number 
of nonlinear effects such as bimodal maser operation^ 
four-wave mixing^ the generation of combs and third 
harmonics,— and a gyrotropic effect.— The interactions in 
the resonator occur between the ions and WG modes, and 
the spatial and frequency distribution of both are funda- 
mental in the nonlinear processes. The ions behave as 
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a number of spin packets with different ESR frequencies 
which overlap to sum to a single broad resonance. The 
greater the homogeneous broadening of the ions, the nar- 
rower the total ESR bandwidth. Thus, it is possible to 
find some spin packets in the crystal having the same 
configuration of frequency and linewidth, but located at 
different sites which make interaction possible with few 
of them. In the section MI Bl the frequency sensitivity 
of the WG modes to a DC magnetic field is developed. 
However, at small magnetic field no description is pre- 
sented. The results of the fit showed in Fig. [5] shows a 
crossing between the two lines different from zero Gauss, 
indicating a shift in the central frequency. For the two 
crystals the shift is 3.7 Gauss, equivalent to 10 MHz shift 
in the central frequency as active internal magnetic field 
is acting. 



IV. ION CONCENTRATION AND SPIN-SPIN 
RELAXATION TIME 

The crystals characterized in the present work are 
HEMEX grade (grown using the heat exchange method) 
single-crystal a-Ak>03. They contain extremely small 
trace concentrations of paramagnetic impurities due to 
the manufacturing process, such as Cr 3+ , Mo 3+ , Ti 3+ 
and Ti 4 +, V+, Mn 2 +, Mn 3 + and Mn 4+ , and Fe+, Fe 2+ , 
and Fe 3+ . The concentration of each ion is approxi- 
mately 1-2 ppm. It has been shown previously through 
measurements of the AC magnetic susceptibility via WG 
modes that the effective concentration of the ion of in- 
terest (Fe 3+ ) is about 10 parts per billion. It has been 
shown that mass conversion of Fc 2+ ions into to Fe 3+ 
can be achieved by annealing the crystal in air, with the 
Fe 3+ concentration increasing by a factor of 10 in crystal 
C2 relative to CI. In parallel to these results it has been 
shown that the spin-spin relaxation time, related to in- 
teractions between nearby spins in the sapphire lattice, is 
about 2 ns - a value smaller than that found in the litera- 
ture of ~10 ns at 4.2 K. This relaxation time is related to 
the ESR linewidth, which is -27 MHz at 4.2 K. Shorter 
spin-spin relaxation times correspond with increased in- 
homogeneous broadening of the ESR linewidth. In this 
section, we characterise this parameter for our crystals as 
a function of the applied magnetic field. The ion concen- 
tration and spin-spin relaxation time can be estimated 
for each crystal from the AC magnetic susceptibility de- 
fined by Eq. The concentration is extracted from the 
DC part of the magnetic susceptibility as follows: 



Xo 



(10) 



where AiVy is the population difference between levels i 
and j, which is determined by a Boltzmann distribution. 
By considering Nij , the total population of the transition 
between i and j, we find 



A Nij = tanh ] 



It is possible to link the frequency v to the corresponding 
magnetic field by a first order approximation of Eq. [3] in 
the limit B < 3000 G: 



v — vb=o i S x B 



v B =a ± S x Bi 



(12) 



Substituting ATVy into the Kramers-Kronig relationships 
gives 



X'(B) = 
X"{B) = 



a 2 xbx S(B-Bij) 
l + {ax S) 2 {B - B^) 2 



a x b 

l + (ax S) 2 (B - B l:j ) 2 



(13) 



where 



2-7TT2 , and 



2h 



fj,oafjNij tanh 



h 



2k B T 



\Vb=o — SBij) 



(14) 
(15) 



(11) 



Note that the parameter a defines the spin-spin relax- 
ation time. From equation [5] it is essential to determine 
the correct filling factor for each WG mode to extract 
the concentration. Due to the number of modes used 
in the fit it is necessary to determine the filling factors 
using a rigorous electomagnctic simulation of the Whis- 
pering Gallery modes.— Because the density of modes is 
high and the simulation process lengthy, the determina- 
tion of the correct modes is often difficult. As a result, 
we extract the product r/N from the real part of the AC 
susceptibility, then normalise the product to the highest 
value obtained from the fit. The result of the fit which 
gives r]N is shown in Fig. [TOJ 

The error bars in the figure are statistical errors from 
the fit results. They are due to the number of points 
around the mode frequency, and to the method of inter- 
action with the mode. The spatial distribution of some 
WG modes does not encompass a large enough quantity 
of ions which makes the detection of the ESR frequency 
difficult. Nonetheless, the results give a good approxima- 
tion to the value of rj with the assumption that a similar 
number of ions interact with all WG modes. The filling 
factors extracted are then used in the imaginary part of 
the AC susceptibility x" t° compute the corrected iV and 
r 2 . The WG modes whose Q-factor degraded by 10% or 
greater when the ESR frequency was coincident with the 
WG mode frequency were selected in the fit (Fig. HU|) . 
For the crystal CI, the concentration is in the range of 10 
to 40 ppb, and for C2 it is in the range 100 to 300 ppb. 
The total concentration should be the same for both Zee- 
man components, as the only difference is the number of 
ions seen by each mode. Therefore the number of ions is 
fixed by the WG mode volumes. For example, the mode 
at 12.04 GHz encloses a volume of about 11.7 x 10~ 6 m 3 , 
which contains a population ofof5.6xl0 15 ions for the 
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crystal CI, and 5.4 x 10 16 ions for the crystal C2. Fig- 
ure [TU] shows the estimated concentrations for the two 
crystals. 
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FIG. 10. Ion concentration, (a) Shows r\ x N in units of 
ion/m 3 , (b) shows the concentration in parts per billion, and 
(c) shows the total concentration calculated for all the modes. 
Only the results for the crystal C2 are shown for clarity. 



The spin-spin relaxation time is deduced from the fit 
in the same way, and the results are summarised in Table 
IIIII A comparison between the real and imaginary parts 
of the AC magnetic susceptibility is made. The results 
are similar to each other and are in close agreement with 
a values found in the literature. 



onance defined by Eq. [16] 



Ai/ = 



1 

7TT 2 



(16) 



The bandwidth of each Zeeman component, and the spin- 
spin relaxation time of the ions are shown in Fig. 1111 
It is clear that the value of 30 MHz is very close to 
the one published by Bogle & Symmons, and Kornienko 
& Prokhorov (see Table HTH . However, some work has 
shown that the ESR bandwidth is concentration depen- 
dent and cannot be accurately measured when highly 
doped (N > 1 ppm)J£ 
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FIG. 11. Spin-spin relaxation time T2 as function of frequency. 
The corresponding ESR bandwidth is also shown. 



TABLE III. Evaluation of the spin-spin relaxation time T2. 
The evaluations for x" were made taking in account the fill- 
ing factor rj. For comparison, the values given by Bogle & 
Symmons^ and Kornienko & Prokhorov^ 3 - are given. 



CI C2 







t 2 (ns) 


At2 (ns) 


r 2 (ns) 


Ar 2 (ns) 




ESRiow 


9.6 


3.6 


12.8 


2.1 


x fit 


ESRhigh 


4.1 


2.7 


7.5 






ESR 


6.9 


4.2 


11.0 


3.2 




ESRiow 


8.1 


3.5 


11.9 


4.8 


x" fit 


ESRhigh 


14.8 


9.5 


11.0 


3.4 




ESR 


10.6 


7.0 


11.4 


3.7 






t 2 (ns) 


At2 (ns) 






Bogle j 


i Symmons 


11.8 


2.2 






Kornienko & Prokhorov 


11.8 


2.4 







Another way of determining t-i is the width of the res- 



V. CONCLUSION 

In conclusion, the physical parameters of Fe 3+ ions in 
high-purity sapphire have been determined at 20 mK by 
measuring the AC susceptibility of the ions using Whis- 
pering Gallery modes. A magnetic field is applied to split 
the ESR of the ion under study into two Zeeman com- 
ponents. The characterisation of the WG modes allowed 
the frequency evolution of each Zeeman component with 
magnetic field to accurately be determined. We measure 
the frequency sensitivities as (2.81 ± 0.022) MHz/Gauss 
for the crystal with high Fc 3+ concentration, and (2.79 ± 
0.014) MHz/Gauss for the crystal with lower concen- 
tration. The zero field splitting has been measured as 
12.039 GHz for both crystals. The difference in terms of 
frequency sensitivity can be explained by an anisotropic 
behaviour of the crystal depending on the concentration 
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of ions involved in the process. However, calculations of 
the g-factor gives results similar to those found in the lit- 
erature. It is expected that the Hamiltonian parameters 
do not change significantly between 4.2 K and 20 mK. 
Indeed, the behaviour of the ESR remains constant for 
temperatures below 4.2 K. To estimate these parameters 
for our crystals it is necessary to know the behaviour 
with magnetic field of the transitions at 19.3 and 31.3 
GHz. Finally, the concentrations estimated arc similar 
to the ones calculated in preceding publication using dif- 
ferent methods, such the AC susceptibility measured by 
saturating the ions at 12.04 GHz, or calculated from a 
Whispering Gallery mode maser oscillator. For the crys- 
tal C2, a concentration of ~200 ppb of Fe 3+ was mea- 
sured, and for CI a concentration of 20 ppb. In previ- 
ous works, the ESR bandwidth of Fe 3+ in sapphire has 



been assumed to be around 27 MHz at 12.04 GHz. Here, 
we measure a value of 34 ± 5.5 MHz for crystal CI at 
12.040 ± 0.0063 GHz, and a value of 29 ± 9 MHz for crys- 
tal C2 at 12.039 ± 0.0038 GHz. Within the bounds of 
experimental error, the bandwidth appears to be inde- 
pendent of the ion concentration. 



ACKNOWLEDGMENTS 

This work was funded by Australian Research Council 
grant numbers DP1092690, DP0986932, CE110001013, 
and FL0992016, and Research Discovery Award grant 
number 12104411. 



* I kar im. benmessai@uwa .edu. aul 

1 P.-Y. Bourgeois, N. Bazin, Y. Kersale, V. Giordano, M. E. 
Tobar, and M. Oxborrow, Appl. Phys. Lett. 87, 224104 
(2005). 

2 K. Benmessai, P.-Y. Bourgeois, Y. Kersale, N. Bazin, M. E. 
Tobar, J. G. Hartnett, M. Oxborrow, and V. Giordano, 
Electron. Lett. 43, 1436 (2005). 

3 P.-Y. Bourgeois, M. Oxborrow, M. E. Tobar, N. Bazin, 
Y. Kersale, and V. Giordano, Int. J. of Mod. Phys. B 20, 
1606 (2006). 

4 K. Benmessai, D. L. Creedon, M. E. Tobar, P.-Y. Bour- 
geois, Y. Kersale, and V. Giordano, Phys. Rev. Lett. 100, 
233901 (2008). 

5 D. L. Creedon, K. Benmessai, M. E. Tobar, J. Hartnett, 
P.-Y. Bourgeois, Y. Kersale, J.-M. Le Floch, and V. Gior- 
dano, IEEE Trans. Ultrason. Ferroelectr. Freq. Control 57, 
282 (2009). 

6 K. Benmessai, M. E. Tobar, N. Bazin, P.-Y. Bourgeois, 
Y. Kersale, and V. Giordano, Phys. Rev. B 79, 174432 
(2009). 

7 K. Benmessai, P.-Y. Bourgeois, M. E. Tobar, N. Bazin, 
Y. Kersale, and V. Giordano, Meas. Sci. Technol. 21, 
025902 (2010). 

8 D. L. Creedon, K. Benmessai, W. P. Bowen, and M. E. 
Tobar, Phys. Rev. Lett. 108, 093902 (2012). 

9 D. L. Creedon, K. Benmessai, and M. E. Tobar, Phys. 
Rev. Lett. 109, 143902 (2012). 

10 G. S. Bogle and H. F. Symmons, Aust. J. Phys. 12, 1 

(1958) . 

11 G. S. Bogle and H. F. Symmons, Proc. Phys. Soc. 73, 531 

(1959) . 

12 H. F. Symmons and G. S. Bogle, Proc. Phys. Soc. 79, 468 
(1962). 

13 L. S. Kornienko and A. M. Prokhorov, Sov. Phys. JETP 
13, 1120 (1961). 

14 K. Benmessai, Maser Cryogenique a Modes de Gallerie , 
Ph.D. thesis, Universite de Franche-Comte (2008). 

15 R. Amsiiss, C. Roller, T. Nobauer, S. Putz, S. Rotter, 
K. Sandner, S. Schneider, M. Schrambock, G. Steinhauser, 
H. Ritsch, J. Schmiedmayer, and J. Majer, Phys. Rev. 
Lett. 107, 060502 (2011). 



T. Wirth, J. Lisenfeld, A. Lukashenko, and A. V. Ustinov, 
Appl. Phys. Lett. 97, 262508 (2010). 

17 P. Bushev, A. K. Feofanov, H. Rotzinger, I. Protopopov, 
J. H. Cole, C. M. Wilson, G. Fischer, A. Lukashenko, and 
A. V. Ustinov, Phys. Rev. B 84, 060501 (2011). 

18 B. Bleaney and R. S. Trenam, Proc. R. Soc. A 223, 15 
(1954). 

19 J. Y. Buzare, G. Silly, J. Klein, G. Scholz, R. Stosser, and 
M. Nofz, J. Phys: Condens Mat 14, 10331 (2002). 

20 S. C. S. C. McGuire, G. P. Lamaze, and E. A. Mackey, 
Trans. Am. Nucl. Soc. 87, 484 (2002). 

21 F. Benabid, M. Notcutt, V. Loriette, L. Ju, and D. G. 
Blair, J. Phys. D: Appl. Phys. 33, 589 (2000). 

22 R. K. Route, M. M. Fejer, A. Alexandrovski, and 
V. Kondilenko, in LIGO LCS meeting , G020374-00-Z 
(2002). 

23 D. L. Creedon, M. E. Tobar, J.-M. Le Floch, Y. Reshitnyk, 
and T. Duty, Phys. Rev. B 82, 104305 (2010). 

24 S. Gekschwind, P. Kisliuk, M. P. Klekin, and D. L. Wood, 
Phys. Rev. 126, 1684 (1962). 

25 R. Borcherts, H. Azarbayejani, S. Karavelas, and G. Sher- 
wood, Microwave Interations at Crystal Defect Centers , 
Report: ORA Project 04275 (College of Engineering, 
Department of Nuclear Engineering, Radiation-Solid-State 
Physics Laboratory, 1963). 

26 T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, 
Phys. Rev. D 77, 053012 (2008). 

27 B. Odom, D. Hanneke, B. D'Urso, and G. Gabrielse, Phys. 
Rev. Lett. 97, 030801 (2006). 

28 R. Zhou, B. Q. Sun, X. Z. Ruan, H. H. Luo, Y. Ji, W. Z. 
Wang, F. Zhang, and J. H. Zhao, J. App. Phys. 103, 
053901 (2008). 

29 L. S. Singer and G. Wagoner, J. Chem. Phys. 37, 1812 
(1962). 

30 J. Stankowski, S. Waplak, and W. Bednarski, Solid State 
Commun. 115, 489 (2000). 

31 P. Pfeffer and W. Zawadzki, J. Appl. Phys. Ill, 083705 
(2012). 

32 J. P. Nibarger, R. Lopusnik, Z. Celinski, and T. J. Silva, 
Appl. Phys. Lett. 83, 93 (2003). 

33 R. M. Eremina, Magnetic Resonance in Solids. Electronic 
Journal 1 (1997). 



10 



A. G. M. A N Luiten and D. G. Blair, J. Appl. Phys. D: 
Appl. Phys. 29, 2082 (1996). 

E. G. Sharoyan, O. S. Torosyan, E. A. Markosyan, and 
V. T. Gabrielyan, Phyica Status Solidi B 65, 773 (1974). 
D. P. Ma, N. Ma, J. Chen, and Q. Lin, Comm. Theor. 
Phys. 33, 167 (2000). 

D. L. Carter and A. Okaya, Phys. Rev. 118, 1485 (1960). 
R. Gallay, J. J. van der Klink, and J. Moser, Phys. Rev. 
B 34, 3060 (1986). 

P. T. Squire and J. W. Orton, Proc. Phys. Soc. 88, 649 
(1966). 



J. Krupka, Derzakowski, M. Tobar, J. Hartnett, and R. G. 

Geyer, Meas. Sci. Technol. 10, 387 (1999). 

M. E. Tobar and J. G. Hartnett, Phys. Rev. D 67, 062001 

(2003). 

J. G. Hartnett, M. E. Tobar, J. M. Le Floch, J. Krupka, 
and P. Y. Bourgeois, Phys. Rev. B 75, 024415 (2007). 
J. le Floch, M. Tobar, D. Cros, and J. Krupka, Appl. Phys. 
Lett. , 58 (2008). 



